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ABSTRACT 


Bubble  formation  during  horizontal  flow  boiling  of  trichlorofluoromethane  (Rll),  1,1- 
dichloro-2,2,2-trifluoroethane  (R123)  and  two  R123/alkylbenzene  lubricant  mixtures  was 
investigated  both  visually  and  calorimetrically.  The  test  fluid  was  pumped  through  the 
inside  of  a roughened,  horizontal,  quartz  tube  which  was  electrically  heated  with  a metal 
strip.  The  refrigerant  entered  the  test  section  with  a quality  slightly  above  the  saturated 
state.  Lxjcally  measured  heat  transfer  coefficients  were  taken  simultaneously  with  high 
speed  motion  picture  images  of  the  boiling  process.  Predictive  equations  from  the 
literature  yielded  acceptable  agreement  with  the  measured  bubble  diameters  and  contact 
angles.  The  addition  of  lubricant  to  the  R123  increased  the  size  of  the  contact  angle  and 
reduced  the  size  of  the  bubble.  The  agreement  between  existing  correlations  and  the 
measured  bubble  frequencies  and  site  densities  was  within  the  uncertainty  of  the 
measurements.  The  addition  of  a small  amount  (0.5%)  of  alkylbenzene  to  R123  increases 
the  number  of  active  nucleation  sites  by  approximately  5 sites/cm^  which  corresponds  to 
a 12%  to  50%  increase  in  the  site  density.  The  increase  in  the  site  density  contributed  to 
the  enhancement  of  the  heat  transfer  coefficient  of  the  R123/0.5%  alkylbenzene  mixture 
over  that  of  the  pure  R123.  Further  increase  in  the  amount  of  alkylbenzene  to  the  R123 
reduces  the  number  of  active  sites  to  below  that  of  pure  R 123  to  approximately  the  value 
for  that  of  Rll.  Consequently,  the  0.5%  lubricant  mass  fraction  mixture  exhibited  a heat 
transfer  coefficient  that  was  larger  than  that  of  the  2%  lubricant  mass  fraction  mixture. 
Correspondingly,  the  lower  heat  transfer  coefficient  of  Rll  as  compared  to  that  of  R123 
was  partially  due  to  the  lower  number  of  active  boiling  sites  for  R1 1.  The  dependency  of 
the  measured  two-phase  heat  transfer  coefficient  on  the  heat  flux  and  Reynolds  number 
was  investigated.  Increases  in  both  the  heat  flux  and  the  Reynolds  number  caused 
increases  in  the  heat  transfer  coefficient.  The  heat  flux  has  a much  larger  effect  on  the 
heat  transfer  coefficient  than  the  Reynolds  number. 
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NOMENCLATURE 


English  symbols 

A cross  sectional  surface  area  of  bubble  (m) 
a major  radius  of  bubble  (m) 

B constant  in  eqn.  19 

b minor  radius  of  bubble  (m) 

Cjj  constant  in  eqn.  19 

Cp  specific  heat  (J/kg  • K) 

D(i  bubble  (departure)  diameter  (m) 

Djj  constant  in  eqn.  19 

Dj  Internal  diameter  of  quartz  tube  (m) 

Dg  equivalent  spherical  diameter  (m) 

E boiling  constant  in  eqn.  15 

f bubble  frequency  (bubbles/s) 

f bubble  frequency  flux  (bubbles/m^  s) 

fps  film  speed  (frames/s) 

g gravitational  acceleration  (m/s^) 

Hj  radial  height  of  image  of  bubble  (m) 

Hq  radial  height  of  bubble  object  (m) 

h2^  two-phase  heat  transfer  coefficient  (W/m^‘K) 

hfg  latent  heat  of  vaporization  (J/kg) 

I constants  in  optical  correction  eqn.  (m^) 

Kp  dimensionless  pressure,  see  eqn.  16 
k thermal  conductivity  (W/m*K) 

m mass  flow  rate  (kg/s) 

Ng  number  of  end-bubbles  for  one  trace 
Nf  number  of  filming  frames 

N_  number  of  zero-bubbles  for  one  trace 
n site  density  (sites/m^^) 

P pressure  (Pa) 

Pe  Peclet  number,  see  eqn.  17 

Pg  critical  pressure  (Pa) 

Pj.  reduced  pressure  = Pg/Pc 

Pg  saturation  pressure  (Pa) 

qg  heat  transfer  due  to  bubble  generation  (W) 

qjj  nucleate  heat  transfer  (bubble  generation  & sensible)  (W) 

q"  heat  flux  (W/m^) 

q"jj  nucleative  heat  flux  (W/m^) 

Re  Reynolds  number  = 4m//xi*'7r*Di 

R^  mean  surface  roughness  (m) 

Rp  Glattungstiefe  or  peak-to-mean  surface  roughness  (m) 

Rq  rms  surface  roughness  (m) 
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r radial  tube  coordinate  with  origin  at  tube  center  (m) 

Sf  standard  deviation  of  measured  bubble  frequency  (s'^) 

s^  standard  deviation  of  measured  jS  (deg.) 

Tg  saturation  temperature  (K) 

temperature  of  the  internal  tube  wall  (K) 
t time  (s) 

Vq  actual  or  object  volume  of  bubble  (m^) 

Vg  volume  of  all  the  end-bubbles  for  one  film-trace  (m^) 
Vj,  volume  of  all  the  zero-bubbles  for  one  film-trace  (m^) 
W mass  fraction  of  lubricant  in  mixture 
X axial  tube  coordinate  (m) 

Greek  symbols 

/3  dynamic  contact  angle  (deg.) 

AT,  (K) 

At  time  difference  between  zero  and  end  time  lines  (s) 

0 angle  between  bubble  major  axis  and  r (deg.) 

/I  dynamic  viscosity  (kg/m  • s) 

V kinematic  viscosity  (m^/s) 

P density  (kg/m^) 

a surface-tension  (N/m) 


Subscripts 

b bubble 

e end 

i image,  inner  tube  surface 

L lubricant 

1 liquid 

m mixture 

o object 

r pure  refrigerant 

s saturation 

V vapor 
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INTRODUCTION 


This  paper  presents  a comparison  of  calorimetric  and  visual  measurements  of  horizontal 
nucleate  flow  boiling  of  four  fluids:  (1)  trichlorofluoromethane  (Rll),  (2)  its  proposed 
replacement,  the  alternative  refrigerant  l,l-dichloro-2,2,2-trifluoroethane  (R123),  (3)  a 
R123/0.5%  weight  alkylbenzene  lubricant  mixture,  and  (4)  a R123/2%  weight 
alkylbenzene  lubricant  mixture;  the  nominal  kinematic  viscosity  of  the  lubricant  was  53 
firn^/s  (280  SUS)  at  313.15  K.  The  calorimetric  aspect  of  this  study  focuses  on  the 
measurement  of  the  local  two-phase  heat  transfer  coefficient  (h2^).  The  visual 
measurements  obtained  from  high-speed  16  mm  film  of  the  boiling  were  taken 
simultaneously  with  the  calorimetric  measurements.  The  bubble  diameters  (D^j),  the 
bubble  frequency  flux  (f),  the  site  density  (n)  and  the  contact  angle  ((3)  were  derived  from 
over  40  high-speed  films  containing  approximately  3000  bubbles. 

Visual  observation  of  the  ebullition  process  is  an  excellent  method  by  which  knowledge 
of  the  fundamental  mechanisms  of  boiling  can  be  obtained.  Currently,  the  heat  transfer 
cannot  be  predicted  a priori  with  a general  mechanistic  model.  In  order  to  establish  the 
foundation  for  this  goal,  fundamental  mechanisms  must  be  observed  and  correlated  for  use 
in  such  a model.  Rohsenow  (1988)  has  brought  to  our  attention  the  need  for  better 
correlations  for  the  bubble  diameters,  the  bubble  frequency  (f),  and  the  site  density.  He 
cited  the  Mikic-Rohsenow  (1969)  model  which  partitions  the  total  heat  flux  between  that 
used  to  generate  the  bubbles  and  that  convected  from  the  surface.  Someday,  fundamental 
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and  generic  mechanistic  boiling  information  will  be  used  to  accurately  design  and  optimize 
refrigeration  equipment.  This  paper  adds  to  the  available  mechanistic  boiling  data  and 
examines  the  applicability  of  several  existing  correlations  to  predict  mechanistic 
information. 

TEST  APPARATUS 


Following  is  a discussion  of  the  test  apparatus  and  the  accuracies  associated  with  the 
individual  measurements.  Unless  otherwise  stated,  all  accuracies  are  estimated  for  a 
99.7%  confidence  interval. 

A schematic  of  the  test  apparatus  is  shown  in  Fig.  1 to  illustrate  the  circulation  of 
refrigerant  through  the  individual  components  of  the  rig.  A hermetic,  oil-free  pump  was 
used  to  deliver  a constant  flow  rate  throughout  the  test  rig.  The  test  fluid  entered  the  fluid 
heater  subcooled  and  exited  in  a two-phase  state  at  near  1%  vapor  quality.  The  electric 
heater  was  used  to  achieve  higher  flow  qualities  entering  the  test  section.  After  exiting 
the  test  section,  the  fluid  entered  a liquid  reservoir  which  served  to  condense  the  vapor 
produced  by  the  quartz  tube  and  the  fluid  heaters  and  to  establish  and  maintain  a steady 
system  pressure.  The  liquid  exiting  the  condenser  was  subcooled  by  approximately  40  K 
and  entered  the  mass  flow  rate  meter.  The  refrigerant  mass  flow  rate  (rh)  was  measured 
with  a Coriolis  flow  meter  with  an  accuracy  of  1-2%  of  the  measured  value.  The  test  loop 
was  completed  by  returning  the  liquid  to  the  pump. 


2 


The  test  section  consisted  of  a 9 mm  internal  diameter  quartz  tube  with  a 3 mm  wall 
thickness.  A thin  (0.25  mm),  3 mm  wide  brass  strip  was  located  horizontally  along  the 
bottom  of  the  tube  with  its  length  aligned  with  the  flow  direction.  A direct  current  was 
passed  through  the  brass  strip  to  produce  a range  of  uniform  heat  fluxes  (q")  from  15  to 
30  kW/m^.  The  quartz  tube  was  enclosed  in  a safety  housing  with  four  flat  windows 
located  at  opposite  poles  for  lighting  and  filming  of  the  boiling.  The  space  between  the 
housing  and  the  quartz  tube  was  evacuated  in  order  to  minimize  the  heat  loss  from  the 
heater  strip  to  the  surroundings  and,  thus,  promote  one-dimensional  conduction  into  the 
tube.  The  all-liquid  Reynolds  numbers  based  on  the  diameter  of  the  quartz  tube  ranged 
from  0.0  to  9,500.  The  accuracy  of  the  Reynolds  number  was  1.4-2. 2%  of  the  calculated 
value.  The  measurements  for  all  fluids  were  obtained  from  the  same  smooth-tube  test  rig 
to  ensure  that  observed  differences  were  due  to  the  characteristics  of  the  fluids. 

The  inside  surface  of  the  quartz  tube  was  roughened  with  a five  micron  polish  to  promote 
boiling.  The  roughness  of  the  tube  was  measured,  using  a stylus  instrument.  The  root 
mean  square  (rms)  roughness  (Rq),  the  Glattungstiefe  (peak- to- mean)  roughness  (Rp),  and 
the  average  roughness  (R^)  of  the  interior  quartz  surface  were  approximately  0.36  ^m,  1 
fxm,  and  0.3  ^m,  respectively.  The  R^  measurement  is  a good  indication  of  the  depth  of 
the  average  cavity.  Close  examination  of  the  roughness  trace  reveals  inclusions  of 
approximately  1 micron  in  width.  Presumably,  the  average  cavity  is  approximately  1 
micron  wide  at  the  mouth  and  approximately  0.3  microns  deep. 
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Care  must  be  taken  in  the  instrumentation  of  a quartz  test  specimen  to  accurately  measure 
the  heat  flux.  The  main  concern  is  determining  what  portion  of  the  energy  generated  in 
the  brass  heater  strip  actually  gets  into  the  tube.  Figure  2 shows  the  instrumentation  of 
the  quartz  tube  with  a heater  composite  which  was  used  to  measure  the  heat  transfer 
coefficient.  The  heater  composite  was  made  by  cementing  two  thermocouples  to  the  center 
of  opposing  sides  of  a plastic  strip  of  known  thermal  conductivity.  The  thermocouple 
nearest  the  tube  was  approximately  0.5  mm  in  diameter  and  flush  with  the  surface  of  the 
plastic.  The  5 mm  x 3 mm  x 76  mm  long  plastic  strip  was  cemented  to  a thin  (0.25  mm) 
brass  strip  of  equal  width  and  length.  The  entire  heater  assembly  was  then  cemented  to 
the  bottom  of  the  quartz  tube.  The  purpose  of  the  assembly  was,  first,  to  direct  most  of 
the  heat  flux  into  the  quartz  tube,  and,  second,  to  provide  for  a method  to  calculate  the 
amount  of  heat  that  enters  the  tube. 

The  heat  flux  incident  to  the  quartz  tube  (q")  and  the  temperature  of  the  inner  tube  wall 
(Twi)  determined  from  both  measurement  and  theory.  For  example,  the  two- 
dimensional  conduction  within  the  plastic  strip  was  modeled  in  order  to  determine  the 
amount  of  energy  transferred  to  the  surroundings  from  the  sides  and  bottom  of  the 
insulating  strip.  The  model  predicted  a linear  heat  loss  with  respect  to  the  temperature 
difference  across  the  plastic.  Consequently,  the  heat  loss  from  the  sides  and  bottom  of  the 
plastic  was  calculated  from  the  measured  temperature  difference  across  the  plastic.  Next, 
the  heat  flux  flowing  into  the  quartz  tube  was  calculated  from  the  difference  between  the 
measured  input  power  to  the  heater  and  the  heat  escaping  from  the  sides  of  the  plastic. 
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The  accuracy  of  the  heat  flux  calculation  was  estimated  to  be  approximately  ±1.7  kW/m^. 
The  estimated  accuracy  of  the  inner  tube  wall  temperature  (T^j)  calculation  was  ± 0.2  K. 
The  T^j  was  extrapolated  from  the  measured  temperature  of  the  outer  brass  strip,  using 
the  calculated  heat  flux  and  a radial  conduction  model. 

The  two-phase  heat  transfer  coefficient  (h2^)  was  calculated  from: 
a" 

uri  s 

where  the  saturation  temperature  (T^)  was  determined  from  the  measured  pressure  to 
within  0.8  K.  The  accuracy  of  the  heat  transfer  coefficient  measurement  for  most  of  the 
data  was  estimated  to  be  ± 10%  of  the  measured  value. 

The  lubricant  mass  fraction  (W)  was  calculated  using  a modification  of  the  procedure  given 
in  ASHRAE  Standard  41.1-84  (ASHRAE  1984).  The  lubricant  mass  fraction  is  the  ratio 
of  mass  of  the  lubricant  to  the  total  mass  of  the  lubricant  and  refrigerant  in  the  mixture. 
The  results  of  three  30-60  g samples  were  averaged.  The  refrigerant  was  evaporated  from 
the  lubricant  by  means  of  a 48  hour  evacuation.  The  mass  fraction  measurement  was 
estimated  to  be  accurate  to  within  0.0005  in  mass  fraction. 
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VISUAL  TECHNIQUE 


A high-speed  16  mm  camera  with  a macro  lens  was  used  to  film  the  boiling  process 
simultaneously  with  the  heat  transfer  measurements.  All  films  were  color  and  taken  at 
approximately  6000  fps.  The  film  speed  was  calculated  with  the  aid  of  timing  marks  on 
the  film  to  within  1.3%  of  the  measured  value.  All  flow  qualities  were  near  1%.  The 
fluid  temperature  entering  the  quartz  tube  that  was  measured  with  a thermocouple  always 
agreed  within  ± 1 K of  the  saturation  temperature  obtained  from  the  corresponding 
measured  pressure.  The  flow  pattern  from  a typical  test  is  presented  in  Fig.  3.  The  liquid 
entering  the  quartz  tube  was  free  of  bubbles  with  the  exception  of  the  large  plug  of  vapor 
that  traveled  along  the  top  of  the  tube. 

The  film  images  were  corrected  for  optical  distortions  induced  by  the  curvature  of  the 
quartz  tube.  Both  the  angle  of  the  major  axis  of  the  bubble  with  respect  to  the  axis  of  the 
tube  (6)  and  its  size  were  altered  by  the  curvature  of  the  tube,  see  Fig.  3.  The  aberrations 
could  have  been  lessened  by  filling  the  cavity  between  the  tube  and  the  protective  housing 
with  a liquid  having  an  index  of  refraction  that  closely  matched  that  of  the  tube.  Doing 
this  would  have  compromised  the  heat  transfer  measurements  since  the  amount  of  heat 
leaving  the  edges  of  the  heater  would  have  been  increased  by  the  presence  of  the  liquid. 
Figure  4 shows  the  optical  correction  curves  that  were  generated  to  correct  the  size  of  the 
bubbles  after  scaling  the  dimensions  from  the  outer  diameter  of  the  tube.  The  ordinate  is 
the  ratio  of  the  radial  height  of  the  image  (H^)  to  that  of  the  actual  object  (Hq)  within  the 
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tube.  The  abscissa  is  the  radial  distance  from  the  center  of  the  tube.  The  real  image  is 
magnified  by  a relatively  constant  value  with  the  exception  of  a small  region  very  close 
to  the  tube  wall  where  it  is  reduced. 


It  was  necessary  to  determine  the  actual  orientation  of  the  bubble  to  calculate  the  actual 
length  of  the  major  radius  (a^,)  of  the  bubble.  The  following  equation  was  developed 
(McClain,  1992)  to  find  the  true  angle  that  the  major  axis  of  the  bubble  made  with  respect 
to  the  tube  axis  (6^): 


0-  = —arccot— — 

° 2 Ij, 

where 


(2) 


Ja  = a^^sin^Q^  + b^^cos^Q^,  (3) 

- a/)sin9^cose^^ 

ia/cos^Qj^  + b/sin^Qj_) 


The  aj  and  the  bj  are  the  minor  and  major  radii  of  the  distorted  bubble.  The  6^  and  the 
are  the  angles  of  the  major  axis  with  respect  to  the  tube  axis  for  the  object  and  the 
image,  respectively.  The  preceding  bubble  dimensions  are  illustrated  in  Fig.  3 where  6-^ 
is  negative  since  it  is  in  the  fourth  quadrant.  The  bubble  image  is  not  rotated  if  Ij,  = 0 
since  the  bubble  cross  section  is  circular. 


7 


The  following  equations  relate  the  actual  major  and  minor  radii  of  the  object  (a^  and  bj 
to  the  radii  of  the  image  and  the  actual  bubble  angle  (McClain,  1992): 


a^b. 

^jgCOS^e^ 

+ J^sine^cosG^  + J^sin^e^ 

/ — • O r\ 

^ . r\  . -r-  Ork 

y/l^sin^d^  - J^sine^cose^  + J^cos^e^ 


One  cautionary  note:  if  the  value  obtained  from  eqn.  7 is  larger  than  that  obtained  from 
eqn.  6 then  eqn.  7 gives  the  value  for  the  major  radius  and  ■7r/2  should  be  subtracted  from 
the  object  angle  which  was  obtained  from  eqn.  2.  Otherwise,  no  corrections  are  necessary 
to  the  above  optical  equations. 


The  volume  and  the  minor  radius  of  the  bubble  were  calculated  from  its  major  radius  and 
the  cross-sectional  area  of  the  actual  bubble  (A^).  The  was  obtained  from  the 
measured  cross-sectional  area  of  the  bubble  image  (A^)  as: 

(8) 


The  major  diameter  was  obtained  by  longitudinally  splitting  the  area  of  the  bubble  in  half 
with  a line  in  the  general  direction  of  the  major  axis  and  measuring  the  distance  from  the 
two  points  of  the  line  that  intersect  the  bubble  edges.  It  was  assumed  that  surface-tension 
and  drag  forces  cause  the  bubble  to  be  symmetric  about  its  major  axis.  Accordingly,  the 
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bubble  volume  was  obtained  by  revolving  the  bubble  about  its  major  axis.  The  ellipsoidal 
volume  was  calculated  from: 


4V 

3713^ 


(9) 


The  accuracy  of  the  volume  calculation  for  two  confidence  intervals  was  estimated  to  be 
within  22%  of  the  measured  value. 

A typical  film-trace  is  depicted  in  Fig.  3.  The  film-traces  were  used  to  calculate  several 
quantities  including  the  amount  of  energy  required  to  generate  the  bubbles  (q^).  Three 
different  types  of  bubbles  were  defined  to  facilitate  the  calculation.  The  first  set  of 
bubbles,  the  zero-bubbles,  are  bubbles  that  are  close  to  the  heated  surface  under  the  zero 
time  line.  The  time  line  is  an  arbitrarily  chosen  reference  line  below  which  all  bubbles 
are  analyzed.  The  time  line  moves  such  that  the  bubbles  that  were  originally  above  the 
zero  time  remain  above  it  and  newly  generated  bubbles  are  below  it.  The  second  set  of 
bubbles,  the  end-bubbles,  are  all  of  the  bubbles  which  are  under  the  time  line  after  it  has 
progressed  from  time-zero  to  some  finite  time  interval  At.  The  difference  between  the 
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volume  of  all  of  the  end-bubbles  and  the  volume  of  all  of  the  zero-bubbles  along  with  the 
time  interval  At  and  the  film  speed  are  used  to  calculate  the  from: 


(V^  - V^)fps 
Nf 


(10) 


The  estimated  accuracy  of  q^,  for  a 68%  confidence  interval,  varied  from  10%  to  80%  of 
the  measured  value,  and  was  typically  within  25  % of  the  measured  value. 

The  third  set  of  bubbles,  the  track-bubbles,  are  bubbles  at  discrete  intervals  of  time  along 
the  time  interval  At.  The  track-bubbles  are  used  to  determine  the  amount  of  heat  the 
bubbles  receive  or  release  after  they  have  departed  from  the  wall.  This  heat  is  used  to 
correct  q^  so  that  only  the  heat  used  to  generate  the  bubbles  (qg)  is  accounted  for.  The 
following  equation  was  developed  to  approximate  the  corrected  nucleative  energy: 


(11) 


where  dq^/dNf  is  the  average  rate  of  change  in  the  energy  received  by  the  bubbles  with 
respect  to  the  number  of  filming  frames.  The  magnitude  of  dq^/dNf  is  equal  to  the  slope 
obtained  from  a linear  regression  of  a plot  of  q^  versus  Nf.  The  correction  was  positive 
approximately  as  many  times  as  it  was  negative.  The  magnitude  of  the  correction  was 
typically  10%  of  q^.  The  sole  purpose  of  the  track-bubbles  was  to  obtain  this  rate  of 
energy  exchange  for  use  in  eqn.  11. 
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CALORIMETRIC  RESULTS 


Pure  Rll  and  Pure  123 

The  two-phase  heat  transfer  coefficients  of  four  different  fluids  were  investigated:  (1)  Rll, 
(2)  R123,  (3)  R123  and  0.55%  alkylbenzene,  and  (4)  R123  and  2%  alkylbenzene.  All 
measurements  were  taken  at  the  same  reduced  pressure  (P^)  of  0.03  and  a flow  quality  of 
near  1%.  Heat  transfer  data  were  taken  for  Reynolds  numbers  varying  from  nearly  0.0 
to  9500  at  three  different  heat  fluxes  (15,  20,  and  30  kW/m^).  Figures  5 through  7 are 
plots  of  the  horizontal  two-phase  heat  transfer  coefficient  versus  the  Reynolds  number  for 
the  three  different  heat  fluxes.  Each  figure  includes  the  data  of  R 123  so  that  a comparison 
between  R123  and  R123/lubricant  mixtures  and  Rll  and  its  replacement  refrigerant  can 
be  readily  made. 

Figure  5 compares  the  heat  transfer  coefficients  of  Rll  to  its  nearly  ozone-safe 
replacement  R123.  The  heat  transfer  coefficient  of  R123  is  larger  than  that  of  R1 1 for  an 
equal  heat  flux  and  Reynolds  number.  At  a zero  Reynolds  number  condition,  the  heat 
transfer  coefficient  for  R123  is  30%  to  42%  greater  than  that  for  Rll.  The  heat  transfer 
advantage  of  R123  steadily  decreases  with  increasing  Reynolds  number.  If  data  for  Rll 
and  R123  are  extrapolated  to  Reynolds  numbers  greater  than  10,000,  one  would  anticipate 
that  the  heat  transfer  coefficient  for  Rll  would  exceed  that  for  R123.  This  is  consistent 
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with  the  measured  data  of  Kedzierski  and  Didion  (1991),  who  observed  that  the  heat 
transfer  coefficient  for  Rll  in  the  nucleate  flow  boiling  regime  is  indeed  higher  than  that 
for  R123  for  Reynolds  numbers  of  18,000  and  24,000. 

Perhaps,  the  mechanism  associated  with  the  elevation  of  the  heat  transfer  coefficient  with 
respect  to  an  increase  in  the  Reynolds  number  for  the  nucleate  flow  boiling  data  presented 
here  is  similar  to  the  mechanism  which  causes  higher  heat  transfer  coefficients  for  tubes 
in  the  upper  rows  of  a smooth-tube  flooded  evaporator.  Stephan  and  Mitrovic  (1981)  have 
found  that  the  convective  enhancement  that  is  present  for  smooth  tube  bundles  is  not 
observed  for  bundles  containing  reentrant  type  tubes.  Unlike  boiling  on  an  enhanced 
surface  where  the  boiling  is  sheltered  from  convection  and  thus  not  dependent  upon  it,  the 
nucleation  on  a smooth  surface  is  exposed  to  the  convection  and  influenced  by  it. 
Consequently,  nucleate  flow  boiling  within  an  enhanced  tube  would  probably  exhibit  less 
of  a Reynolds  number  dependence  than  that  shown  here. 

The  heat  transfer  coefficient  increases  nearly  linearly  with  respect  to  Reynolds  number  for 
both  fluids.  The  data  for  each  heat  flux  for  each  fluid  are  nearly  parallel  to  one  another, 
indicating  that  the  rate  of  increase  in  h^p  with  respect  to  Re  is  independent  of  the  incident 
heat  flux.  Also,  notice  that  the  Reynolds  number  has  more  of  an  effect  on  the  heat 
transfer  of  R1 1 than  it  does  on  R123.  For  example,  the  two-phase  heat  transfer  coefficient 
for  Rll  increases  by  an  average  of  37%  for  an  increase  in  Re  from  0 to  9500.  For  the 
same  increase  in  the  Reynolds  number,  the  R123  heat  transfer  increases  by  only  10%. 
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For  both  Rll  and  R123,  the  heat  flux  has  a much  greater  effect  on  the  heat  transfer 
coefficient  than  the  Reynolds  number.  For  example,  a 33  % increase  in  the  heat  flux  from 
15  to  20  kW/m^  corresponds  to  an  average  increase  in  the  heat  transfer  coefficient  of  22% 
for  both  Rll  and  R123.  Likewise,  a 50%  increase  in  the  heat  flux  from  20  to  30  kW/m^ 
corresponds  to  an  average  increase  in  the  heat  transfer  coefficient  of  56%  for  Rll  and 
46%  for  R123.  Clearly,  this  boiling  phenomenon  is  heat  flux  driven,  which  is  what  one 
would  expect  given  it  is  for  a constant  heat  flux  boundary  condition. 

R123  and  Alkvlbenzene  Mixture 

It  was  assumed  that  the  change  in  the  boiling  point  of  R 123  due  to  the  addition  of  the 
lubricant  was  negligible.  This  assumption  is  justified  by  two  arguments.  First,  Stephan 
(1963)  has  reported  that  the  boiling  point  of  R12  increases  by  less  than  0.3  K with  an 
addition  of  10%  oil  by  weight.  The  lubricant  mass  fractions  measured  here  are  2%  and 
0.5%  which  should  cause  a comparatively  small  change  in  the  boiling  point.  Second,  no 
change  in  the  agreement  between  the  measured  saturation  temperature  and  the  saturation 
temperature  obtained  from  the  measured  pressure  was  observed  with  the  addition  of 
alkylbenzene  to  the  R123. 

The  properties  of  the  refrigerant/lubricant  mixture  were  evaluated  at  the  saturation 
temperature  of  the  pure  refrigerant.  The  liquid  density  of  the  lubricant  (pl)  the 
kinematic  viscosity  of  the  lubricant  (vi)  were  fitted  to  limited  data  (Chevron,  1992)  as  a 


13 


function  of  temperature  in  kelvins.  The  liquid  density  of  the  lubricant  in  kilograms  per 
cubic  meters  was  fitted  to  the  form: 

= 919.69  + 0.1994r^  - O.OOIBIST/  (12) 

The  kinematic  viscosity  of  the  alkylbenzene  lubricant  in  square  meters  per  second  was 
fitted  to  the  form: 

= 5 .6054xl0~^^exp  ^308.70 


The  following  ideal  mixing  rule  was  used  to  approximate  the  dynamic  viscosity  of  the 
refrigercuit/lubricant  mixture: 

= explfvTnp^  + (l-h^)lnp^  (14) 

The  surface-tension  between  the  liquid  lubricant  and  air/lubricant  vapor  (aj^)  was 
calculated  from  the  rise  height  in  a 1.2  mm  bore  capillary  tube.  The  surface-tension  was 
found  to  be  approximately  0.026  N/m  at  296.4  K.  The  surface-tension  of  R123  was 
obtained  from  Chae  et  al.  (1990).  The  surface-tension  of  Rll  was  obtained  from  the 
International  Institute  of  Refrigeration  (1982). 

Figure  6 compares  the  measured  two-phase  heat  transfer  coefficient  of  pure  R123  to  a 
mixture  of  R123  and  0.55  weight  percent  of  alkylbenzene.  The  heat  transfer  coefficient 
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is  plotted  versus  the  liquid  Reynolds  number.  Two  aspects  of  the  comparative  heat 
transfer  are  evident  from  the  figure.  First,  the  addition  of  0.5%  alkylbenzene  to  pure 
R123  enhances  the  heat  transfer  for  all  heat  transfer  conditions.  Overall,  the  heat  transfer 
coefficient  of  the  R123/0.5%  mixture,  averaged  for  all  Reynolds  numbers  and  heat  fluxes, 
is  16%  greater  than  the  overall  averaged  heat  transfer  coefficient  for  pure  R123.  Second, 
the  addition  of  the  lubricant  to  R123  has  caused  the  heat  transfer  to  be  more  responsive 
to  an  increase  in  the  Reynolds  number.  For  example,  the  heat  transfer  coefficient  for  the 
R123/0.5%  alkylbenzene  mixture  increases  roughly  by  44%  and  11%  for  an  increase  in 
the  Reynolds  number  from  zero  to  9,500  for  the  two  lower  heat  flux  conditions  and  the 
highest  heat  flux  condition,  respectively.  As  a result,  the  enhancement  of  the  R123/0.5% 
alkylbenzene  mixture  heat  transfer  relative  to  that  of  pure  R123  increases  with  Reynolds 
number. 

Figure  7 compares  the  measured  two-phase  heat  transfer  coefficient  of  pure  R123  to  a 
mixture  of  R123/2.  % weight  percent  of  alkylbenzene.  The  heat  transfer  coefficient  of  the 
R123/2%  alkylbenzene  mixture  can  be  examined  in  terms  of:  (1)  its  relative  magnitude 
compared  to  the  pure  component,  and  (2)  the  influence  of  the  Reynolds  number  upon  it. 
Examination  of  the  first  point  above  reveals  that  the  heat  transfer  coefficient  of  the 
R123/2%  alkylbenzene  mixture,  averaged  over  all  heat  transfer  conditions,  is  1 1 % greater 
than  that  for  pure  R123.  This  is  approximately  two-thirds  of  the  enhancement  associated 
with  the  0.5%  mass  fraction  fluid.  A study  of  the  second  point  above  indicates  that  the 
rate  of  increase  for  the  R123/2%  alkylbenzene  mixture  is  greater  than  that  for  pure  R123 
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and  it  is  a also  a function  of  the  heat  flux.  For  example,  an  increase  in  the  Reynolds 
number  from  zero  to  9,500  causes  a 75%,  26%  and  a 23%  increase  in  the  heat  transfer 
coefficient  for  the  15  kW/m^,  20  kW/m^,  and  30  kW/m^  heat  flux  conditions, 
respectively.  Averaged  over  all  heat  transfer  conditions,  the  percent  increase  in  the  heat 
transfer  coefficient  with  respect  to  increased  Reynolds  number  for  the  2 % mass  fraction 
mixture  is  more  than  twice  as  large  as  it  is  for  the  0.5%  mixture. 

Perhaps  the  increase  in  the  rate  of  increase  in  the  heat  transfer  coefficient  with  respect  to 
the  Reynolds  number  is  related  to  the  increase  in  the  liquid  thermal  conductivity  and  the 
liquid  Prandtl  number.  The  thermal  conductivity  of  the  2%  alkylbenzene  mixture  is  0.2% 
greater  than  that  of  the  pure  fluid.  The  Prandtl  number  of  the  2%  mixture  is  12%  greater 
than  that  of  the  pure  fluid. 

Correlation  of  Heat  Transfer  Coefficient 

It  is  difficult  to  find  a heat  transfer  correlation  that  applies  strictly  to  horizontal  nucleate 
flow  boiling.  One  candidate  correlation  form  might  be  a superposition  model. 
Historically,  these  correlations  are  generated  from  a statistical  fit  of  many  flow  boiling 
data  sources  from  many  different  metal  tubes.  Typically,  the  effect  of  surface  roughness 
is  not  included  in  the  correlation  of  a superposition  model.  Yet,  as  Corty  and  Foust 
(1955)  have  shown,  pool  boiling  strongly  depends  upon  the  surface  finish.  Therefore, 
rather  than  comparing  the  present  data  with  a superposition  model  for  flow  boiling,  the 
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present  data  for  a zero  Reynolds  number  are  compared  to  three  different  pool  boiling 
correlations  which  account  for  surface  roughness  effects. 

Agreement  between  the  measured  heat  transfer  coefficients  and  two  of  the  correlations  was 
not  acceptable.  On  average,  the  heat  transfer  coefficient  for  Rll  obtained  from  Cooper’s 
equation  (1984)  was  overpredicted  by  approximately  300%.  Likewise,  his  equation 
overpredicted  the  heat  transfer  coefficient  for  R123  by  approximately  200%.  Similarly, 
poor  agreement  was  achieved  with  Stephan’s  and  Abdelsalam’s  (1980)  pool  boiling 
equation  for  refrigerants  with  the  correction  for  surface  roughness. 


Good  agreement  was  achieved  between  the  measured  heat  transfer  coefficient  at  Re  = 0 
and  the  correlation  with  the  Borishanskii-Minchenko  equation  (Borishanskii  et  al.,  1969): 


N 


9^(Pi  - Pv) 


where  Kp  is  the  dimensionless  pressure: 


= 


a (pj  - pj 


(16) 
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and  the  Peclet  number  (Pe)  is: 


Pe 


q"  pj  c 


pj 


hfg  Pv  ^^Pi  ■ Pv) 


(17) 


Equation  15  and  the  measured  heat  transfer  coefficients  for  Rll  and  R123  at  Re  = 0 
agreed  to  within  ±12%  when  the  dimensionless  boiling  constant  E was  set  equal  to  2.5 
X 10*^.  Rice  and  Calus  (1972)  correlated  the  pool  boiling  data  from  an  electroplated 
chromium  surface  of  Cichelli  and  Bonilla  (1945)  to  the  Borishanskii-Minchenko  equation 
with  the  constant  E = 3.92  x 10*^.  Although  roughness  values  for  this  surface  were  not 
given,  a typical  R^  roughness  for  an  electroplated  surface  is  approximately  0.8  /xm 
(Baumeister,  1978).  If  the  constant  for  Cichelli’s  and  Bonilla’s  data  were  to  be  used  to 
predict  the  present  Rll  and  R123  data,  the  data  would  be  overpredicted  by  an  average  of 
60%.  However,  this  value  of  E may  be  inappropriate  for  the  quartz  tube  since  the 
roughness  of  the  chromium  surface  is  greater  than  that  of  the  quartz  tube.  Consequently, 
the  value  of  E = 2.5  x 10*^  for  the  quartz  tube  seems  to  be  consistent  with  Cichelli’s  and 
Bonilla’s  data  since  one  would  expect  that  the  value  of  E for  the  quartz  tube  would  be 
smaller  than  that  for  the  chromium  surface. 

VISUAL  RESULTS 


Visual  data  with  an  acceptable  statistical  certainty  is  difficult  to  obtain  for  two  reasons:  (1) 
the  cost  of  frequent  repetition  of  visual  measurements  is  prohibitive,  and  (2)  the  time  for 
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which  data  is  taken  is  of  the  order  of  milliseconds.  If  the  boiling  process  is  cyclic,  the 
problem  of  a limited  data  sample  is  worsened. 

Vapor  production  during  nucleate  flow  boiling  can  vary  with  respect  to  time  for  what  is 
viewed  calorimetrically  as  steady  state  heat  transfer.  Strenge  et  al.  (1961)  and  Schultz  and 
Cole  (1978)  have  noted  the  fluctuating  behavior  of  nucleate  boiling.  Mesler  and  Banchero 
(1958)  noticed  that  during  the  boiling  of  organic  liquids,  the  surface  temperature 
fluctuations  decreased  for  higher  pressures  and  higher  heat  fluxes.  The  temperature 
fluctuations  were  not  noticeable  for  pressures  above  atmospheric  pressure. 

All  of  the  data  presented  here  are  for  pressures  slightly  above  atmospheric.  The  hope  is 
that  the  amplitude  of  the  fluctuations  is  small  since  the  data  are  taken  for  a small  time 
interval  (0.025  seconds).  However,  the  amplitude  of  the  fluctuation  can  be  measured  only 
with  repeated  measurements.  If  the  amplitude  of  the  fluctuation  is  small,  then  averaging 
a large  amount  of  data  is  not  as  important.  It  is  likely  that  an  adequate  average  for  the 
bubble  diameter  and  the  contact  angle  can  be  obtained  since  these  have  sufficient  numbers 
for  averaging.  However,  the  site  density  and  the  bubble  frequency  are  more  strongly 
dependent  upon  the  variations  of  the  ebullition  cycle.  Consequently,  as  Kedzierski  and 
Didion  (1990)  have  found,  the  fluctuations  of  the  bubble  frequency  and  the  site  density  are 
greater  than  those  of  the  bubble  diameter  and  contact  angle.  Therefore,  a cautious  view 
should  be  taken  toward  the  bubble  frequency  and  the  site  density  data  for  a particular  heat 
flux  and  Reynolds  number.  For  most  cases,  the  data  have  been  gathered  from  one  film 
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with  a duration  of  approximately  one  second.  In  other  words,  the  frequency  and  site 
density  data  for  a particular  heat  transfer  condition  represent  an  event  which  may  not 
necessarily  be  the  most  probable  event,  i.e.,  the  average  event. 

Foaming  occurred  for  the  zero  Reynolds  number  condition  for  all  of  the 
refrigerant/lubricant  mixtures.  Foaming  caused  nearly  a 50%  increase  in  the  bubble 
diameter  and  approximately  a 50%  reduction  in  the  bubble  flux.  Consequently,  only  the 
non-zero  Reynolds  number  data  for  the  bubble  flux  or  the  bubble  diameter  were  included 
in  the  calculation  of  the  overall  averages  so  that  a fair  comparison  between  each  fluid 
could  be  made. 

Contact  Angle 

The  dynamic  and  the  static  contact  angles  of  all  of  the  fluids  were  measured.  All  test 
fluids  fully  wetted  a clean  glass  surface  when  placed  on  the  surface  as  liquid  droplets. 
Wetting  implies  that  the  static  contact  angle  is  zero.  The  dynamic  contact  angles  were 
obtained  from  the  film-traces.  Figure  3 defines  the  contact  angle  (/?)  as  it  was  measured 
here. 

Figure  8 shows  the  measured  contact  angle  for  the  four  fluids  versus  Reynolds  number. 
Approximately  ten  contact  angles  were  average  from  each  film-trace.  A linear  line  fitted 
through  the  data  would  indicate  that  the  contact  angle  increases  slightly  with  respect  to  the 
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Reynolds  number.  This  makes  intuitive  sense  since  if  a bubble  is  not  abruptly  stripped 
from  the  wall  by  drag  forces  it  will  grow  more  slowly  and  tend  to  flatten  out  resulting  in 
a smaller  contact  angle.  However,  an  average  value  for  the  contact  angle  may  prove  more 
convenient  since  the  influence  of  the  Reynolds  number  on  the  contact  angle  is  relatively 
small. 

The  present  measured  average  contact  angles  confirm  the  approximate  contact  angle  for 
refrigerants  of  35°  cited  by  Stephan  and  Abdelsalam  (1980).  Row  1 of  table  1 presents 
the  contact  angles  for  each  fluid  averaged  over  all  heat  transfer  and  flow  conditions.  For 
example,  the  average  measured  contact  angle  for  pure  Rll  was  31°  which  is  only  12% 
greater  than  the  that  given  by  Stephan  and  Abdelsalam  (1980).  The  average  measured 
contact  angle  for  pure  R123,  36°,  agrees  closely  with  35°.  The  addition  of  alkylbenzene 
to  R123  resulted  in  a larger  contact  angle.  Specifically,  the  average  contact  angle  for 
R123  and  0.55%  and  2%  alkylbenzene  were  40°  and  41°,  respectively. 

Row  2 of  table  1 presents  the  estimated  standard  deviation  of  the  contact  angle 
measurement  (s^).  A portion  of  the  deviation  of  the  measurement  may  be  due  to  the 
variation  of  the  contact  angle  while  the  bubble  grows. 
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Bubble  Diameter 


The  bubble  diameter  appeared  to  be  independent  of  the  Reynolds  number  and  the  heat 
flux.  The  cyclic  characteristics  of  boiling  should  not  have  a significant  effect  on  the  size 
of  the  bubble.  Surface-tension  and  the  contact  angle  should  be  the  primary  factors  in 
determining  the  bubble  diameter.  Therefore,  given  the  above  two  conditions,  an  average 
of  all  the  bubbles  for  a particular  fluid  should  result  in  a fair  representation  of  the  average 
bubble  diameter  for  that  fluid  at  a given  thermodynamic  state. 

An  equivalent  spherical  bubble  diameter  (D^)  was  calculated  from  the  average  volume  of 
the  zero-bubbles.  The  average  diameter  of  the  zero-bubbles  should  closely  represent  the 
departure  diameter.  The  accuracy  of  the  bubble  diameter  measurement  for  two  confidence 
intervals  was  estimated  to  be  within  7%  of  the  calculated  value.  Row  3 of  table  1 presents 
the  average  equivalent  spherical  diameters  for  the  four  fluids  for  non-zero  Reynolds 
numbers.  Row  4 of  table  1 shows  the  standard  deviation  from  the  average  bubble 
diameter  for  each  fluid. 

Fritz  (1935)  utilized  the  work  of  Bashforth  and  Adams  (1883)  to  calculate  the  maximum 
volume  of  a spherical  vapor  bubble  as  a function  of  Laplace’s  constant  and  the  contact 
angle  /3  (in  degrees)  with  respect  to  the  solid  surface.  Fritz’s  solution  is  presented  in 
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graphical  and  tabular  form.  The  tabulated  solution  can  be  fitted  to  the  following  familiar 
equation  for  the  departure  diameter  attributed  to  Fritz: 


0.0204P 


N 


a 

£r(Pi  - Pv) 


(18) 


The  Fritz  equation  does  a remarkable  job  of  predicting  the  data  for  such  a simple  equation. 
When  the  measured  contact  angle  and  the  fluid  properties  are  substituted  into  eqn.  18,  the 
agreement  between  the  equation  and  the  measured  bubble  diameter  for  R123  is  within  1 %. 
Equation  18  underpredicts  the  bubble  diameter  for  Rll  by  approximately  16%.  However, 
the  Fritz  equation  does  not  predict  the  effect  of  the  addition  of  lubricant  to  the  refrigerant. 
Nevertheless,  it  predicts  the  measured  bubble  diameter  for  the  R123/0.55%  alkylbenzene 
and  the  R123/2%  alkylbenzene  mixtures  within  14%  and  35%,  respectively. 


Site  Density 

The  determination  of  the  site  density  was  subjective.  It  was  difficult  to  determine  exactly 
how  many  sites  were  active  since  some  sites  were  very  close  to  each  other.  A best  guess 
was  made  concerning  the  average  number  of  locations  from  which  bubbles  originated. 
The  surface  area  was  calculated  from  the  product  of  the  measured  width  of  the  heater  and 
the  scaled  length  of  the  analyzed  region  obtained  from  the  film-trace. 
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Figure  9 shows  the  heat  flux  dependency  on  the  nucleation  site  density  for  all  four  fluids. 
The  symbols  represent  the  measured  data.  The  solid  lines  are  best  fit  regressions  of  the 
data  to  the  form: 

q"  = ^ (19) 

where  the  constant  Djj  was  set  to  zero  for  all  of  the  fluids  with  the  exception  of  the 
R123/2%  alkylbenzene  mixture.  A linear  fit  was  best  for  this  data.  The  B exponents  for 
the  other  R123  fluids  were  also  close  to  1,  namely,  0.91  for  R123  and  1.24  for  R123  and 

0. 5%  alkylbenzene.  The  exponent  for  the  site  density  for  the  Ril  was  0.46  which  was 
smaller  than  that  for  the  R123  fluids.  The  values  for  B are  summarized  in  row  5 of  table 

1. 

The  B exponent  for  the  fluids  tested  here  does  not  vary  significantly  from  that  reported  in 
the  literature.  Hsu  and  Graham  (1976)  remark  that,  although  the  constant  Cjj  varies 
greatly  depending  upon  the  fluid  and  the  surface,  the  exponent  has  a relatively  small  range 
of  variance  between  about  0.3  and  0.8.  Kirby  and  Westwater  (1965)  have  measured  an 
exponent  for  the  nucleation  site  density  on  a glass  surface  and  found  it  to  be  approximately 
equal  to  0.73. 

The  addition  of  a small  amount  of  alkylbenzene  (0.5  %)  to  R123  increases  the  number  of 
active  nucleation  sites  by  approximately  5 sites/cm^.  Further  increase  in  the  amount  of 
alkylbenzene  to  the  R123  reduces  the  number  of  active  sites  to  approximately  the  value  for 
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that  of  Rll.  For  the  same  heat  flux,  Rll  has  approximately  10  less  active  sites/cm^  than 
those  of  R123.  The  sixth  row  of  table  1 summarizes  the  site  density  measurements. 

Bubble  Flux 

The  bubble  frequency  flux  f"  is  the  total  number  of  bubbles  emitted  from  a surface  per 
unit  area  per  unit  time.  The  bubble  frequency  flux  was  calculated  as: 

f"  = (20) 

Pv  ^fg 

The  accuracy  of  the  bubble  flux  calculation  was  estimated  to  be  approximately  33%  for 
two  confidence  intervals. 

It  appears  that  the  bubble  flux  is  approximately  the  same  for  all  fluids.  Row  7 of  table  1 
gives  the  bubble  frequency  flux  averaged  over  all  non-zero  Reynolds  number  heat  transfer 
conditions  for  each  of  the  fluids.  The  estimated  standard  deviation  of  the  bubble  flux 
measurement  (s^),  given  in  row  8,  is  nearly  greater  than  the  measurement.  Therefore, 
the  bubble  fluxes  for  each  fluid  cannot  be  given  a strict  relative  ranking.  The  bubble  flux 
is  enhanced  for  small  lubricant  mass  fractions  (0.5%)  with  caveat  that  the  standard 
deviation  of  this  measurement  is  large. 
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The  average  bubble  frequency  for  a single  bubble  is  obtained  by  multiplying  the  frequency 
flux  by  the  site  density.  Rows  9 and  10  of  table  1 gives  the  measured  frequency  and  the 
estimated  standard  deviation  of  the  measurement,  respectively.  Based  on  the  work  of 
Jakob  (1949),  Zuber  (1959)  developed  an  equation  for  the  frequency  of  a single  bubble: 


f = 


0 . 59  ^ g 9"  (Pj  - Py) 


(21) 


Predicted  values  of  f from  eqn,  21  lie  within  the  fluctuations  of  the  measured  f.  The 
equation  underpredicts  the  measured  average  bubble  frequency  (f)  by  approximately  100% 
for  all  of  the  test  fluids.  For  example,  eqn.  21  predicts  the  bubble  frequency  to  be  74 
bub/s  for  Rll;  the  measured  value  averaged  for  all  conditions  is  178  bub/s.  The  standard 
deviation  of  the  individual  measurements  from  this  mean  value  is  approximately  182  bub/s. 
Similar  results  were  obtained  for  the  other  fluids  which  are  presented  in  row  9 and  10  of 
table  1. 


MECHANISTIC  INTERPRETATION  OF  ^2^ 

The  visual  data  were  taken  with  a view  to  gaining  some  physical  understanding  of  the 
calorimetric  data.  The  model  of  Mikic  and  Rohsenow  (1969)  can  be  used  as  a tool  for 
making  the  physical  connection  between  the  visual  and  the  calorimetric  data.  Their  model 
of  the  nucleate  heat  flux  (q"jj)  is: 
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= 2 41  n (22) 

This  model  is  based  on  the  assumption  that  the  latent  heat  is  not  important  in  determining 
the  heat  transfer  due  to  nucleation.  Instead,  the  governing  mechanism  is  the  transport  of 
the  superheated  liquid  layer  from  the  wall  by  the  motion  of  the  bubble.  Consequently, 
eqn.  22  is  based  on  a transient  conduction  analysis  on  an  unit  area  which  is  twice  the 
bubble  diameter.  Lorenz  et  al.  (1974)  have  validated  this  model  for  water  and  methanol 
pool  boiling. 

The  heat  transfer  coefficient  obtained  from  equation  22  consists  of  the  multiplication  of 
a fluid  property  part  (2{kjpjCpi7r}°'^)  and  a bubble  mechanism  part  (f''^  n).  The  second 

to  the  last  row  in  table  1 shows  that  the  fluid  property  terms  for  all  of  the  fluids  are  within 
1%  of  each  other.  Accordingly,  the  model  suggests  that  the  observed  differences  in  the 
heat  transfer  coefficients  for  each  fluid  cannot  be  attributed  to  the  differences  in  fluid 
transport  properties.  Consequently,  the  differences  in  the  heat  transfer  coefficients  from 
fluid  to  fluid  must  be  due  to  differences  in  the  bubble  parameters. 

The  last  row  of  table  1 gives  the  magnitude  of  the  bubble  mechanism  term  calculated  using 
the  average  values  for  f,  D^j  and  n which  are  given  in  table  1 . Comparison  of  the  relative 
magnitudes  of  the  bubble  mechanism  term  for  each  fluid  shows  that  this  term  is  larger  for 
the  R123/0.5%  alkylbenzene  mixture  primarily  due  to  the  larger  number  of  active 
nucleation  sites  associated  with  that  fluid.  The  bubble  diameter  and  the  bubble  frequency 
have  not  changed  much  with  the  addition  of  a small  amount  of  lubricant.  Consequently, 
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part  of  the  reason  that  the  R123/0.5%  alkylbenzene  heat  transfer  is  greater  than  that  of 
pure  R123  is  that  the  addition  of  lubricant  has  caused  more  sites  to  become  active  in 
generating  bubbles.  On  the  other  hand,  further  increase  in  the  lubricant  mass  fraction  has 
decreased  the  bubble  size  and  the  site  density  which  has  contributed  to  a reduction  in  the 
heat  transfer  as  compared  to  the  0.5%  lubricant  mass  fraction.  Likewise,  the  lower  heat 
transfer  coefficient  of  pure  Rll  as  compared  to  that  of  pure  R123  is  primarily  due  to  a 
lower  number  of  active  sites. 

CONCLUSIONS 


Bubble  formation  during  horizontal  flow  boiling  of  Rll,  R123  and  a R123/alkylbenzene 
mixture  was  investigated  both  visually  and  calorimetrically.  The  two-phase  heat  transfer 
coefficient  of  four  different  fluids  were  investigated;  (1)  Rll,  (2)  R123,  (3)  R123  and 
0.55%  alkylbenzene,  and  (4)  R123  and  2%  alkylbenzene.  Good  agreement  was  achieved 
between  the  measured  heat  transfer  coefficients  for  Rll  and  R123  and  the  Borishanskii- 
Minchenko  equation. 

The  relative  magnitudes  of  the  nucleate  flow  boiling  heat  transfer  coefficient  were 
compared  to  that  of  R123.  For  Reynolds  numbers  below  9,500,  the  heat  transfer 
coefficient  of  R123  was  on  average  22%  greater  than  that  of  Rll  for  equal  heat  flux  and 
Reynolds  number.  The  addition  of  alkylbenzene  to  pure  R123  enhances  the  heat  transfer 
relative  to  that  of  pure  R123  for  all  heat  transfer  conditions. 
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For  both  Rll  and  R123,  an  increase  in  either  the  heat  flux  or  the  Reynolds  number 
increases  the  heat  transfer  coefficient.  The  heat  flux  has  a much  larger  effect  on  the  heat 
transfer  coefficient  than  the  Reynolds  number.  An  increase  in  the  heat  flux  caused  nearly 
an  equal  percentage  increase  in  the  heat  transfer  coefficient.  The  two-phase  heat  transfer 
coefficient  for  Rll  increases  significantly  more  than  that  of  R123  for  an  increase  in 
Reynolds  number.  The  addition  of  lubricant  to  the  R123  caused  the  heat  transfer  to  be 
more  sensitive  to  an  increase  in  the  Reynolds  number. 

The  average  measured  contact  angle  for  pure  R123  agrees  closely  with  that  reported  in  the 
literature:  35°  as  compared  to  36°.  The  addition  of  alkylbenzene  to  R123  resulted  in  a 
larger  contact  angle.  Namely,  the  average  contact  angle  for  R123  and  0.55%  and  2% 
alkylbenzene  was  40°  and  41°,  respectively. 

Given  the  simplicity  of  the  Fritz  equation  (eqn.  18),  it  does  a remarkable  job  of  predicting 
the  measured  bubble  diameters:  within  1%  and  16%  of  the  measured  values  for  the  pure 
R123  and  Rll,  respectively.  The  addition  of  alkylbenzene  caused  a reduction  in  the 
average  size  of  the  bubble,  which  is  not  predicted  by  the  Fritz  equation.  Even  so,  it 
predicts  the  measured  bubble  diameter  for  the  R123/0.55  % alkylbenzene  and  the  R123/2% 
alkylbenzene  mixtures  within  14%  and  35%,  respectively. 

The  addition  of  a small  amount  (0.5%)  of  alkylbenzene  to  R123  increases  the  number  of 
active  nucleation  sites  by  approximately  5 sites/cm^.  Further  increase  in  the  amount  of 
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alkylbenzene  to  the  R123  reduces  the  number  of  active  sites  to  approximately  the  value  for 
that  of  Rll.  For  the  same  heat  flux,  Rll  has  approximately  10  less  active  sites/cm^  than 
that  of  R 123. 

Zuber’s  equation  (eqn.  21)  predicts  the  measured  average  bubble  frequency  (f)  to  within 
the  measurement  uncertainty  for  all  of  the  test  fluids.  The  total  flux  of  bubbles  from  the 
surface  is  approximately  the  same  for  all  of  the  fluids. 

Part  of  the  reason  that  the  R123/0.5%  alkylbenzene  heat  transfer  is  greater  than  that  of 
pure  R123  is  that  the  addition  of  lubricant  has  caused  more  sites  to  become  active  in 
generating  bubbles.  The  bubble  diameter  and  the  bubble  frequency  has  not  changed  much 
with  the  addition  of  a small  amount  of  lubricant.  On  the  other  hand,  further  increase  in 
the  lubricant  mass  fraction  has  decreased  the  bubble  size  and  the  site  density  which  has 
contributed  to  a reduction  in  the  heat  transfer.  Primarily  a lower  number  of  active  sites 
has  caused  the  pure  Rll  heat  transfer  to  be  lower  than  that  of  the  pure  R123. 
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Table  1 - Visual  Data  Summary 


Measured 

Data 

Rll 

R123 

R123/0.5% 

alkylbenzene 

R 123/2% 
alkylbenzene 

j8  (deg.) 

31 

36 

40 

41 

Sq  (deg.) 

6 

5 

6 

9 

Ds  (mm) 

0.83 

0.76 

0.74 

0.65 

Sds 

0.31 

0.28 

0.28 

0.19 

B 

0.46 

0.91 

1.24 

1 

^ 1 21.7  kW/m2(Cf^*^) 

20 

30 

40 

20 

f (s'^m-2) 

380  X 10^ 

390  X 10^ 

490  X 10^ 

370  X 10^ 

ap(s’^m’^) 

290  X 10^ 

190  X 10^ 

250  X 10^ 

270  X 10^ 

f (s*^) 

178 

144 

146 

251 

Sf  (s'^) 

182 

54 

75 

164 

2{kiP|Cp,x}“'^ 

1185 

1172 

1174 

1179 

(s‘*^) 

1.84 

2.08 

2.65 

1.34 
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© Temperature  measurement 


© Pressure  measurement 
M Valve 


Fig.  1 Schematic  of  test  apparatus 


Fluid,  Ts 


Fig.  2 Schematic  of  heater  composite 


Fig.  3 Trace  method  and  flow  pattern 
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of  tube  tube  wall 


Fig.  4 Optical  correction  ratio  for  radial  height 
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Fig.  5 Two-phase  heat  transfer  coefficient  measurements  for  pure  Rll  and  pure  R123  at 
Pr  = 0.03 
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ha*  (kW  / m^K 


Fig.  6 Two-phase  heat  transfer  coefficient  measurements  for  R123/0.5%  alkylbenzene 
mixture  and  pure  R123  at  Pj.  = 0.03 
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Fig.  7 Two-phase  heat  transfer  coefficient  measurements  for  R123/2%  alkylbenzene 
mixture  and  pure  R123  at  Pj.  = 0.03 


43 


P(deg.) 


60 


50 


40 


30 


20  - 


10  - 


8 


8 


o 

A 

A 

s 


□ Pure  R123 
O Pure  R1 1 

O R123/2%  Alkylbenzene 
A R1 23 /.5%  Alkylbenzene 

_J ^ ^ ^ 


0 


3 4 

Re  X 1 0^ 


Fig,  8 Measured  contact  angle  for  Rll,  R123  and  R123  and  alkylbenzene 
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Fig.  9 Heat  flux  dependency  on  the  nucleation  site  density  for  Rll,  R123  and  R123  and 
alkylbenzene 
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